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LONG  TERM  GOALS 

Time  and  space  dependent  radiance  distributions  at  the  sea  surface  are  a  function  of  the  shape  of  the 
incident  distribution  on  the  surface,  modification  by  the  sea  surface  itself  from  topography  and 
transmission  characteristics,  and  alteration  by  the  Inherent  Optical  Properties  (IOPs)  of  the  surface 
ocean.  Our  long  term  goal  is  understanding  this  last  controlling  factor.  With  a  knowledge  of  the  IOPs, 
radiance  fields  can  be  directly  computed  from  the  incident  field  using  the  equation  of  radiative  transfer, 
now  embedded  in  commercially  available  code  (e.g.,  Hydrolight). 

With  the  state  of  current  technology  and  methodologies,  the  primary  obstacles  in  understanding 
subsurface  IOPs  and  their  high-frequency  dynamics  are  a  lack  of  1)  volume  scattering  instrumentation, 
2)  comprehensive  inversion  models  linking  the  IOPs  with  the  ambient  particle  fields  including  bubbles 
(models  which  in  many  cases  will  require  input  dependent  on  1),  and  3)  suitably  stable,  non-intrusive 
platforms  to  sample  the  subsurface  ocean.  The  first  two  challenges  are  addressed  in  this  project. 

OBJECTIVES 

There  are  two  broad  objectives  for  this  project: 

1)  To  develop  an  in-situ  volume  scattering  function  device  measuring  volume  scattering  from  10° 
to  170°  at  10°  intervals  and  sampling  rates  of  1  s'1  or  better  to  sample  the  VSF  in  near-surface 
waters;  and 
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2)  To  develop  and  refine  IOP  inversion  models  to  resolve  particle  field  characteristics  on  small 
spatial  (cm’s)  and  temporal  (<1  s)  scales  in  near-surface  waters. 

APPROACH 

Our  design  for  the  VSF  device  is  illustrated  in  Figure  1.  The  device  is  called  MASCOT  (Multi-Angle 
Scattering  Optical  Tool).  The  source  beam  is  a  30  mW  658  nm  laser  diode  expanded  with  a  Gallilean 
2X  beam  expander  to  an  approximately  3  mm  X  8  mm  elliptical  shape.  A  wedge  depolarizer  is  used  to 
provide  the  unpolarized  light  needed  for  VSF  determinations.  Seventeen  independent  silicon  diode 
detectors  spaced  in  a  semicircle  10  cm  around  the  sample  volume  measure  volume  scattering  from  10° 
to  170°  at  10°  intervals.  The  total  pathlength  for  all  scattering  measurements  is  20  cm.  Independent 
detectors  allow  resolution  of  the  VSF  without  any  moving  parts  and  time-consuming  scanning. 
Additionally,  each  detector  can  be  optimized  for  its  specific  dynamic  range.  Detector  field-of-views 
(FOVs)  range  from  0.8°  to  5°  for  the  different  detectors,  with  the  narrowest  FOVs  associated  with  the 
detectors  measuring  scattering  in  the  forward  direction.  Using  proprietary  electronics,  a  20  Hz 
sampling  rate  for  all  channels  has  been  achieved  while  maintaining  a  worst  case  signahnoise  of  300: 1 . 
Relatively  fast  sampling  rates  are  important  in  resolving  VSFs  in  the  highly  dynamic  ocean  subsurface. 

Polarized  VSFs  have  now  been  successfully  collected  with  the  addition  of  filter  mount  placed  in  front 
of  the  source  beam.  A  linear  polarizer  is  used  to  obtain  scattering  from  a  vertically  and  horizontally 
polarized  source  (in  terms  of  the  Mueller  scattering  matrix  elements,  (Si  l+S12)/2  and  (Si  1-S  12)/2  are 
measured,  so  that  SI  1  and  S12  may  be  derived).  Adding  polarized  scattering  increases  the  amount  of 
information  on  particle  characteristics  we  are  collecting,  and  is  expected  to  improve  our  ability  to 
discriminate  different  particle  types  (both  the  number  of  subpopulations  and  the  accuracy  of  individual 
determinations).  The  degree  of  linear  polarization  (S12/S1 1)  is  most  dependent  on  the  degree  of 
sphericity,  particle  size,  and  refractive  index. 

For  VSF  inversion  modeling,  we  are  extending  the  capabilities  of  existing  models  (Twardowski  et  al. 
2001;  Twardowski  and  Zaneveld,  2004;  Zhang  et  al.  2005)  by  incorporating  input  from  new  VSF 
measurements  and  by  adding  bubble  particle  populations  (clean  and  coated)  in  the  models.  Candidate 
phase  functions  for  particle  subpopulations  are  fit  to  measured  VSFs  using  a  least-squares 
minimization  matrix  inversion  procedure.  These  phase  functions  can  be  obtained  theoretically  using 
Mie  theory,  DDA,  or  IGOM  techniques,  or  experimentally  in  controlled  laboratory  conditions.  George 
Kattawar’s  group  has  recently  provided  married  DDA-IGOM  phase  function  determinations  for 
asymmetric  polyhedrons  we  are  now  using  in  the  inversions. 

The  MASCOT  has  now  been  extensively  deployed  concurrently  with  the  commercially  available  near¬ 
forward  VSF  device  LISST  (Sequoia  Inc.)  in  order  to  capture  the  VSF  with  good  resolution  from  ~0. 1 
degrees  to  170  degrees.  Our  deployment  package  also  contains  CTD,  AC9  (or  ACS),  and  various  ECO 
scattering  sensors,  all  integrated  in  a  vertical  profiling  system  (Figure  2). 

Deployments  for  RaDyO  have  taken  place  in  collaboration  with  other  RaDyO  investigators  off  Scripps 
Pier  in  January  2008,  in  Santa  Barbara  Channel  in  September  2008,  and  off  Hawaii  in  September  2009. 
For  these  RaDyO  deployments,  we  have  also  intermittently  integrated  two  additional  sensors:  1)  a  high 
sampling  rate  fish-eye  lens  radiometer  currently  being  developed  by  Marlon  Lewis  and  Scott  McLean 
of  Satlantic,  and  2)  a  bubble  acoustic  resonator  developed  by  Svein  Vagle,  David  Farmer,  and  Helen 
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Czerski.  For  the  Scripps  Pier  experiment,  a  mobile  deployment  platform  was  configured  so  that 
measurements  could  be  made  along  the  length  of  the  pier,  with  the  expectation  that  higher 
concentrations  of  bubbles  would  be  observed  nearing  the  surfzone.  For  subsequent  field  work  in  SBC 
and  Hawaii,  we  are  deploying  our  sensor  package  off  the  R/V  Kilo  Moana. 


Figure  1.  Oblique  view  illustration  of  the  MASCOT.  The  VSF  is  resolved  from  10  to  170  degrees  in 
10  degree  intervals.  Detectors  are  wedge  shaped  and  arranged  in  a  semi-circle  on  an  aluminum 
frame  to  minimize  reflections  and  perturbation  of  the  water  sample  in  the  remote  volume  ( center  of 
semi-circle).  The  source  assembly  includes  a  30  m  W  658  nm  laser  diode,  reference  detector,  beam 
expander,  and  wedge  depolarizer.  Wiring  from  all  the  detector  modules  and  the  source  module 

feeds  to  a  data  handling  unit. 


WORK  COMPLETED 

•  An  IOP  package  including  the  MASCOT  and  LISST  was  deployed  off  Scripps  Pier,  in  the 
Santa  Barbara  Channel,  and  off  Hawaii  during  collaborative  RADYO  exercises. 

•  The  MASCOT  was  adapted  to  additionally  measure  linear  polarized  scattering  elements  to 
derive  degree  of  linear  polarization  (DOLP)  as  a  function  of  angle. 

•  An  acoustic  bubble  resonator  (Vagle  and  Farmer)  and  fish-eye  lens  camera  (Lewis)  were 
successfully  integrated  on  the  MASCOT  IOP  package  and  deployed  for  all  exercises. 

•  The  MASCOT  IOP  package  was  tested  off  the  New  York  Bight  in  May  and  July  2009  in 
association  with  other  funded  work. 

•  A  revised  calibration  protocol  was  developed  and  implemented,  using  clean  techniques  and 
accurately  accounting  for  pure  water  absorption  along  the  optical  path. 

•  A  revised  calibration  and  correction  protocol  was  developed  for  the  LISST  to  account  for 
optical  geometry,  ambient  light  contamination,  turbulent  scattering,  and  vignetting. 

•  A  second  MASCOT  prototype  with  a  larger  dynamic  range  suitable  for  high  turbidity 
environments  (such  as  the  surf  zone)  was  tested  in  the  lab  and  deployed  in  the  NY  Bight  and 
the  surf  zone  at  Duck,  NC. 


3 


•  Inversion  of  MASCOT  VSFs  collected  off  Scripps  Pier,  SBC,  and  Hawaii  into  component 
VSFs  associated  with  subpopulations  of  the  bulk  particle  field  including  bubbles  has  been 
refined  using  a  least-squares  minimization  matrix  inversion  model,  including  derivation  of 
particle  concentration  and  size  distribution  for  each  subpopulation.  The  model  has  been 
adapted  to  also  invert  for  nonspherical  asymmetric  polyhedron  particles  more  analogous  to 
natural  suspended  mineral  populations. 

•  Based  on  results  from  this  project,  4  papers  have  been  published  by  our  group,  1  additional 
paper  has  been  submitted,  2  additional  papers  are  close  to  submission,  and  2  additional  papers 
are  in  preparation  for  the  RaDyO  JGR  Special  Issue. 

RESULTS 

Fig.  2  shows  the  MASCOT  device  in  a  custom  cage  with  additional  optical  sensors  ready  for 
deployment.  Orientation  in  the  horizontal  plane  minimizes  any  shear  perturbation  of  the  sampled 
water  parcel  during  upcasts. 

We  have  implemented  the  VSF  inversions  now  with  phase  functions  computed  for  nonspherical 
particles  (asymmetric  polyhedrons).  These  have  been  assimilated  into  our  phase  function  library  along 
with  Mie  theory  derived  phase  functions  for  coated  bubbles.  Phase  functions  for  asymmetric 
polyhedron  monodispersions  are  shown  in  Fig.  3. 

Inversion  results  from  the  Scripp’s  Pier  surf  zone  using  combinations  of  asymmetric  polyhedron 
subpopulations  (mineral-mimicking),  coated  spherical  particles  (bubbles)  and  Very  Small  Particles 
(VSPs)  are  shown  in  Fig.  4.  The  VSPs  are  quasi-Rayleigh  scatterers  with  sizes  near  the  wavelength  of 
the  source  beam;  for  these  particles,  composition  information  is  difficult  to  glean  from  scattering 
measurements  alone.  Inversion  results  are  substantially  more  stable  using  the  asymmetric  polyhedron 
phase  functions  as  opposed  to  previous  phase  functions  where  nonspherical  particles  were 
approximated  with  spheres.  Particles  dominated  by  sediments  were  present  over  the  first  hour  of  the 
time  series,  with  bubbles  becoming  much  more  important  as  the  surf  zone  approached  our  sampling 
area  with  the  ebbing  tide.  Fig.  5  shows  the  inverted  subpopulation  size  distributions  during  the  time 
period  of  localized  wave  breaking.  This  is  the  first  time  the  theoretically  expected  mid-angle 
enhancement  in  scattering  due  bubbles  has  been  observed  in-situ.  The  20  Hz  sampling  rate  for  all  VSF 
channels  allows  sufficient  resolution  of  the  bubble  plumes  as  they  rapidly  evolve  in  time  and  space.  A 
very  low  background  scattering  condition  (particle  attenuation  less  than  0.3  m"1  at  532  nm)  enhanced 
our  dynamic  range  in  resolving  intense  scattering  associated  with  episodic  bubble  and  sediment  plume 
generation. 

Inversion  results  are  encouraging  and  are  consistent  with  anecdotal  evidence  (what  we  saw)  with 
respect  to  patches  of  sediment  and/or  bubbles.  Fig.  6  shows  the  inversion  results  for  the  2  primary 
sediment  subpopulations,  and  Fig.  7  shows  results  for  the  primary  large  bubble  subpopulation.  Fig.  8 
evaluates  the  use  of  P(70)/(3(  120)  as  a  proxy  for  large  bubble  subpopulation  abundances.  The  positive 
correlation  is  encouraging,  although  more  work  is  required  to  see  if  using  specific  bins  in  the  larger 
size  classes  may  improve  the  fit.  Fig.  9  shows  acoustics  data  collected  with  the  Svein/Farmer  bubble 
resonator  concurrently  with  the  large  bubble  subpopulation  from  the  inversion  for  the  last  40  minutes 
of  the  Fig.  6  time  series.  Very  good  agreement  between  patterns  in  optical  scattering  from  bubbles 
(positive  displacements  of  the  p(70)/p(120)  relationship  relative  to  the  background)  and  acoustic 
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attenuation  at  multiple  frequencies  was  observed.  Fig.  10  shows  the  size  distributions  from  the 
inversion  results  for  one  of  the  measured  VSFs,  with  an  aggregate  distribution  consistent  with  the 
familiar  Junge-type  distribution  in  the  optically  relevant  size  range  (~0.1  to  -40  pm).  Stable  model 
solutions  (only  possible  to  assess  with  a  very  large  number  of  VSFs)  consistent  with  anecdotal 
observations  of  particle  composition,  coupled  with  validation  based  on  acoustic  measurements  of 
bubble  distributions  and  fractional  and  aggregate  size  distributions  agreeing  with  theory  and 
observations,  provides  strong  evidence  that  our  inversion  results  are  environmentally  meaningful. 

We  also  able  to  leverage  other  ongoing  work  to  deploy  the  MASCOT  device  along  with  the 
Farmcr/Vaglc/Czcrski  acoustic  resonator  in  East  Sound,  WA  in  May  2010.  This  exercise  helped 
evaluate  the  importance  of  surfactants  in  stabilizing  bubble  populations  through  coating  of  the  bubbles. 
For  the  experiment,  a  collaborating  boat  passed  off  our  bow  at  full  throttle  creating  a  linear  propeller 
wash.  Our  boat  then  towed  the  MASCOT  and  resonator  back  and  forth  through  the  propeller  wash 
along  with  a  holographic  camera  from  Joe  Katz’s  group  (JHU).  Jim  Chumside  also  participated  in  the 
experiment,  flying  lidar  missions  over  the  propeller  wash.  Substantial,  repeatable  signatures  from 
bubbles  were  observed  in  the  wash  by  the  optical  and  acoustical  techniques,  with  persistence  for  over  1 
hour  (Fig.  11).  Similar  measurements  from  the  day  before  in  waters  with  low  surfactants  revealed  that 
propeller  wash  bubbles  vanished  after  10-15  minutes.  Surfactants  are  thus  very  important  in 
stabilizing  near-surface  bubble  populations. 

Preliminary  inversion  results  from  Hawaii  for  a  time  series  where  the  acoustic  bubble  resonator  was 
also  deployed  with  the  MASCOT  are  shown  in  Fig.  12.  Interestingly,  there  is  no  evidence  of  a 
persistent  small  bubble  population  in  this  data  set,  as  was  observed  in  the  inversion  results  from 
Scripp’s  Pier.  During  bubble  injections  and  subsequent  seconds  of  bubble  plume  evolution,  the  peaks 
in  bubble  size  distributions  were  typically  observed  between  1-2  pm. 

General  observations  from  our  VSF  inversions  and  validation  efforts  can  be  summarized  as: 

•  Patterns  in  bubble  and  mineral  contributions  indicate  consistently  stable  solutions  when 
nonspherical  particles  are  considered  (high  number  of  VSFs  essential  to  assess  this) 

•  Patterns  agree  qualitatively  with  anecdotal  observations 

•  Concurrent  acoustic  measurements  corroborate  large  bubble  inversion  results 

•  PSDs  from  inversion  results  follow  commonly  observed  Junge-type  distribution  reasonably 
well 

•  When  modeling  all  particles  as  spheres,  inversion  is  less  stable,  but  most  mineral  and  bubble 
subpopulation  trends  are  accounted  for 

IMPACT/APPLICATIONS 

Progress  and  results  represent  important  steps  toward  the  development  and  vetting  of  a  multi-angle,  in¬ 
water  VSF  device.  Knowledge  of  the  Inherent  Optical  Properties  including  the  VSF  can  be  used  to 
predict  and  optimize  the  performance  of  a  host  of  Naval  operations  that  rely  on  divers,  cameras,  laser 
imaging  systems,  and  active  and  passive  remote  sensing  systems.  These  include  mine 
countermeasures,  harbor  security  operations,  debris  field  mapping,  anti-submarine  warfare,  and  search 
and  salvage  operations. 
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Figure  2.  MASCOT  VSF  device  mounted  with  other  optical  sensors  in  a  custom  cage 
designed  to  sample  the  subsurface  domain. 
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Figure  3.  Phase  functions  of  asymmetric  polyhedrons  computed  from  a  combination 
of  the  Discrete  Dipole  Approximation  (DDA)  and  Improved  Geometric  Optics  Model 
(IGOM)  to  span  a  wide  range  of  particle  sizes.  Calculations  were  carried  out  in  George 
Kattawar’s  lab  in  collaboration  with  Yu  You  and  Li  Bi.  Asymmetric  polyhedrons  were 
generated  by  slightly  adjusting  the  axes  of  a  symmetric  rectangular  polyhedron  to 
mimic  the  crystalline  structure  of  mineral  quartz. 
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Figure  4.  Time  series  from  VSF  measurements  off  Scripp ’s  Pier.  Over  time,  the  surf  zone  moved 
from  a  shoreward  location  to  the  vicinity  of  our  sensing  package  with  an  ebbing  tide.  The  VSF 
resolved  at  120  degrees  is  plotted  in  the  upper  panel.  In  the  first  hour,  4  scattering  events,  marked 
as  “ wave  sets ”  in  the  panel,  were  detected.  Each  involved  a  remnant  sediment  plume  from 
shoreward  surf  zone  resuspension  that  became  entrained  in  a  rip  current  running  along  the  pier. 
Approximately  18  min  separation  between  the  plumes  corresponded  with  primary  wave  sets  hitting 
the  surf  zone.  After  the  first  hour,  enhanced  scattering  is  due  to  localized  resuspension  and  bubble 
entrainment  from  breaking  waves.  The  ratio  of  scattering  at  70  deg  to  scattering  at  120  deg  is 
shown  in  the  middle  panel.  Bubble  particles  exhibit  enhanced  scattering  between  60  and  80  deg,  so 
this  ratio  is  a  sensitive  indicator  of  bubbles.  Within  the  first  hour,  bubbles  are  only  detected  by  this 
proxy  ratio  in  a  random  local  breaking  event  at  36  min.  Afterward,  substantial  peaking  of  this  ratio 
coincides  with  localized  breaking  and  bubble  injection.  The  bottom  panel  shows  background 
particle  population  results  from  our  VSF  inversion  for  aggregate  bubble  populations  (blue), 
aggregate  mineral  subpopulations  (cyan)  and  “  Very  Small  Particles  ”  that  cannot  be  distinguished  in 
terms  of  composition.  Mineral  particles  are  observed  to  dominate  the  background  population,  with 
some  evidence  for  a  persistent  small  bubble  population  centered  at  about  0.2  pm. 
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Figure  5.  Same  as  Fig.  4  except  bottom  panel  shows  inverted  size  distributions  within 
a  scattering  peak  dominated  by  bubbles  from  local  injection.  The  distribution  for  minerals 
exhibits  a  broadening  into  larger  size  classes.  The  primary  difference  in  these  populations 
relative  to  the  background,  however,  is  a  very  substantial  bubble  population  that  dominates 

the  larger  particle  sizes  (>  20  pm). 
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Figure  6.  3  panel  plot  showing  the  ratio  of  VSF  at  70  deg  to  VSF  at  120  deg,  and  the  abundances 
for  the  two  dominant  sediment  subpopulations  over  time.  The  bottom  panel  is  the  larger 
distribution.  Substantial  peaks  are  observed  during  the  4  sets  within  the  first  hour  due  to  sediment 
plumes  entrained  in  a  rip  current.  Afterward,  substantial  abundances  are  observed  in  localized 

resuspension  from  wave  breaking. 
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Figure  7.  Plot  showing  the  ratio  of  VSF  at  70  deg  to  VSF  at  120  deg  (upper  plot) 
and  the  dominant  large  bubble  subpopulation  (lower  plot).  In  the  first  hour,  larges  bubbles  are  only 
observed  during  a  random  local  breaking  event  at  36  min.  Substantial  abundances 
of  large  bubbles  are  observed  afterward  in  associated  with  localized  wave  breaking 

and  bubble  injection. 
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N  (m"  )  in  large  bubble  subpopulation 


Figure  8.  Large  bubble  subpopulation  abundances  from  the  VSF  inversion  plotted  against  the 
simple  ratio  of  the  VSF  at  70  deg  to  VSF  at  120  deg.  A  linear  correlation  is  observed  with  an  r2  of 
0.67,  indicating  that  the  ratio  may  be  used  as  a  simple  proxy  to  approximate  large  bubble 

distributions  over  time. 
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Figure  9.  Inversion  results  from  the  large  bubble  subpopulation  (upper  panel)  plotted  with  relative 
acoustic  attenuation  from  a  bubble  resonator  (lower  panel).  Patterns  were  very  consistent,  more 
consistent  than  observed  when  comparing  the  acoustic  attenuation  with  positive  responses  in 
(l(70)/ji(l20)  front  the  MASCOT.  The  time  series  shown  here  is  a  subset  of  the  data  shown  in  Figs. 

4-7,  including  only  wave  sets  4  and  5.  Note  that  the  turbidity  maximum  observed  around  58-60 
minutes  in  Fig.  6  is  not  associated  with  any  response  from  the  acoustics  or  the  inverted  large  bubble 
subpopulation.  This  is  consistent  with  the  previous  interpretation  that  this  turbidity  plume  is 
dominated  by  sediments.  A  particularly  close  agreement  was  observed  between  the  inversion  results 
and  acoustic  attenuation  at  the  frequency  specific  for  143  um  diameter  bubbles.  Note  that  even 
though  the  two  sensors  were  on  the  same  instrument  package,  sample  volumes  were  separated  by 

about  60  to  70  cm. 
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Figure  10.  Size  distributions  of  dominant  subpopulations  with  relative  refractive  indices  of  0. 75 
(bubbles)  and  1.16  (quartz-like  minerals),  and  Very  Small  Particles  (VSPs).  The  aggregate  size 
distribution  (in  black)  is  modeled  well  with  a  hyperbolic  differential  Junge-type  model  with  slope  3. 7 
in  the  range  that  can  be  resolved  by  sizing  devices  such  as  a  Coulter  Counter.  Many  previous 
studies  looking  at  oceanic  size  distributions  have  observed  that  the  Junge  distribution  is  often  a 
reasonable  model,  which  serves  as  an  element  of  validation  for  our  inversion  results. 
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Figure  11.  Optical  and  acoustical  evidence  of  bubble  persistence  behind  a  boat’s  propeller 
wash  for  over  an  hour  in  waters  front  East  Sound,  WA  that  exhibited  high  surfactant 
content.  Previous  data  collected  in  waters  with  low  surfactant  showed  bubble  populations 
rapidly  vanishing  after  10-15  min.  Data  show  the  importance  of  surfactant  coatings  in 
stabilizing  near-surface  bubble  populations. 
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Figure  12.  Time  series  of  scattering  at  120  deg  (upper  panel),  ratio  of  scattering  at 
70  deg  to  scattering  at  120  deg  (middle  panel),  and  PSD  inversion  results  for  minute  26 
(lower  panel)  for  Hawaii  data  collected  09/2009.  Phase  function  library  for  the  inversion 
only  included  a  background  phase  function  and  phase  functions  for  bubble  subpopulations 
computed  from  Mie  theory.  Interestingly,  there  is  no  evidence  of  a  persistent  small 
bubble  population  in  this  data  set,  as  was  observed  in  the  inversion  results  from  Scripp ’s 
Pier.  During  bubble  injections  and  subsequent  seconds  of  bubble  plume  evolution,  the 
peaks  in  bubble  size  distributions  were  typically  observed  between  1-2  pm. 

TRANSITIONS 

We  expect  that  our  efforts  in  developing  an  in- water  VSF  device  and  associated  inversion  techniques 
to  better  understand  particle  dynamics  in  natural  waters  will  lead  to  transition  as  operational  tools  for 
the  fleet  and  the  oceanographic  research  community  in  the  future. 
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RELATED  PROJECTS 


This  effort  is  related  to  several  ongoing  efforts  by  the  PI  to  develop  optical  sensors  and  associated 
biogeochemical  inversion  techniques  to  improve  our  understanding  of  the  oceanic  environment. 

Current  ongoing  related  projects  include: 

•  investigating  the  dynamics  of  scattering  by  subsurface  bubble  populations  and  other  particles  in 
the  S.  Ocean  (NASA,  Twardowski;  project  lead  PI  H.  Dierssen); 

•  investigating  the  underlying  controls  of  biological  camouflage  responses  in  dynamic  underwater 
optical  environments  (MURI  collaboration,  Twardowski;  URI  Pis  J.  Sullivan  and  B.  Seibel; 
project  lead  PI  M.  Cummings); 

•  developing  improved  remote  sensing  water  quality  algorithms  for  coastal  waters  (NASA, 
Twardowski;  project  lead  PI  Z-P.  Lee); 

•  developing  compact,  low  power  sensing  tools  for  ocean  observing  platforms  (ONR  SBIR, 
Twardowski); 

•  developing  a  microscopic  holographic  camera  for  optically  relevant  particles  (NOPP,  project 
lead  PI  J.  Sullivan); 

•  developing  optical  prediction  models  for  the  surfzone  (CEROS,  Twardowski;  project  lead  PI  G. 
Chang);  and 

•  developing  a  surfzone  drifter  measuring  optical  attenuation  and  scattering  (ONR  SBIR, 
Twardowski). 
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